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We present results for 1 — 2 pion scattering length with the Wilson fermions in the quenched approximation. 
The finite size method presented by Liischer is employed, and calculations are carried out at /3 = 5.9, 6.1, and 
6.3. In the continuum limit, we obtain a result in reasonable agreement with the experimental value. 



1. Introduction 

Lattice calculations of scattering lengths of the 
two-pion system is an important step for under- 
standing of strong interactions beyond the hadron 
mass spectrum. For the / = process, which is 
difficult due to the presence of disconnected con- 
tributions, only one group carried out the calcu- 
lation [01 . For the 1 = 2 process, on the other 
hand, a number of calculations has been car- 
ried out with the Staggered and the Wil- 
son fermions ||l|,||. These calculations reported 
results in agreement with the prediction of cur- 
rent algebra or lowest-order chiral perturbation 
theory(CHPT) [|. It is known, however, that 
this prediction differs from the experimental value 
over 1.4(7, and that the higher order effects of 
CHPT are small @. 

The past lattice calculations were made on 
coarse lattices with small sizes, and the contin- 
uum extrapolation was not taken. In this arti- 
cle we report on our high statistics calculation 
of the 7 = 2 pion scattering length aiming to 
improve on these points. This work is carried 
out in quenched lattice QCD employing the stan- 
dard plaquette action for gluons with the Wilson 



fermions. The number of configurations (and lat- 
tice size) are 187(163 x 64), 120(24^ x 64), and 
100(323 X 80) for /3 = 5.9, 6.1, and 6.3. The pion 
mass covers the range 450 — 900MeV. 

2. Method 

The energy eigenvalue of a two-pion system in 
a finite periodic box jg shifted by finite-size 
effect. Liischer presented a relation between the 
energy shift AE and the S'-wave scattering length 
ao given by 



AE- 



T + A-T'^ + B-T^ + 0{T^) , (1) 



'presented by N. Ishizuka 
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where T = ao/(7rL). Since T takes a small 
value, typically ^ — lO^'^, in our simulation we 
can neglect the higher order terms 0{T^). The 
constants A and B are geometrical values A = 
-8.9136 • • • and B = 62.9205 • • •. 

The energy shift AE can be obtained from the 
ratio R{t) = G{t)/D(t), where 

G{t) = {7T+it)7T+it)W-{h)W-{t2)) 

D{t) = {n+{t)W-{h)) {n+{t)W-{t2)) • (2) 

In order to enhance signals we use wall sources 
(denoted by W~) and fix gauge configurations to 
the Coulomb gauge. The two wall sources are 
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Figure 1. R{t) = G{t)/D{t). 

placed at different time slices ti and ^2 to avoid 
contaminations from Fierz-rearranged terms in 
the two-pion state that would occur for the choice 
ti — t2. In this work we set ^2 = ti + 1 and 
ti = 8, 10, 13 for (3 = 5.9, 6.1, 6.3. Quark prop- 
agators are solved with the Dirichlet boundary 
condition imposed in the time direction and the 
periodic boundary condition in the space direc- 
tions. In region t >> ti,t2 the ratio behaves as 
R{t) - Z{1 ~AE-{t~ ti) + 0{t^)). 

As an example, the ratio R{t) at /3 = 6.3 and 
K = 0.1513 corresponding to = 433(4)MeV 
is plotted in Fig. |l|. The signal is very clear and 
Z ^ \. This means the overlap of our wall sources 
with the two-pion state is sufficiently large. In 
general there are higher order terms 0{t^) in 
but we cannot resolve them in Fig. 0. Making a 
linear fitting in the range i = 27 — 62, we obtain 
(aA£;) = 5.97(60) x 10"^. Solving equation (§) 
with this value we obtain T = -1.73(15) x 10"^, 
which corresponds to Oq = — 0.525(45)(l/GeV). 

3. Result 

In Fig. H we compare our new results for 
ao/a™^"'" with those of old calculations, where 
is the prediction of current algebra : 
= — m7r/(167r/^). For the decay constant 
we use the value at finite to^ at finite lat- 
tice spacing referred in each study. This ratio 
has been commonly employed to make a com- 
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Figure 2. Comparison of our results with the old 
calculations by a^/a^^^^ . 



parison of current algebra and lattice calculations 
with different quark actions and parameters. The 
open symbols refer to results with the Staggered 
fermions and filled ones are those of the Wilson 
fermions. The legends give /3, the spatial lattice 
size L, and collaborations of the studies ( K : Ku- 
ramashi et.al. §, SGK : Sharpe et.al. ||, GPS : 
Gupta et. al. Q , Ours : our calculation ) . We also 
plot the experimental value at ra^ = 140MeV. 
We find that our results are inconsistent with 
old results, especially with those of the Staggered 
fermions. 

A possible cause of the discrepancy is the sys- 
tematic error of determination of /^r needed to 
calculate Cq^^""". In Fig. ^ we compare our re- 
sults with old calculations in terms of ao/m^. 
The same symbols as those in Fig. ^ are used. 
The lattice results including ours are almost con- 
sistent with each other. Also they appear to be 
in more agreement with the experiment than with 
the prediction of current algebra. 

We note that the calculation of being de- 
termined by the amplitude of correlation function 
of pion and axial vector current, is quite difficult. 
Various systematic errors may well enter in their 
determinations. Further the mass dependence of 
is not small and oq /oq ^^"^ is very sensitive to 
it. For these reasons we analyze ao/m^ below. 
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Figure 3. Comparison of our results with the old 
calculations by ao/mTr- 

From chiral symmetry oq /tott behaves as 

ao/m^ = A+B-ml+C-ml\og{ml/A^)+0{m'^) .(3) 

For the Wilson fermions we should consider an- 
other term oc l/rn^ that arises from breaking of 
chiral symmetry. Golterman and Bernard also 
proposed the same term pointing out that it can 
appear from quenching effect Q]. However, these 
effects are very small in our simulation as we do 
not observe a rapid variation of af)/m„ expected 
from such a term in Fig. ^. Further the chiral 
logarithm term, log(m^/A^), is also small. 

In Fig. ^ our results for ao/wTr in the chiral 
limit obtained by a linear fitting in are plot- 
ted, together with the experimental value and the 
prediction of current algebra at m^r = 140MeV. 
We observe a very clear linear dependence in the 
lattice spacing a. By a linear extrapolation, we 
then obtain 



ao/m^ = -1.91(25) (l/GeV^) 
aom„ = -0.0374(49) , 



(4) 



in the continuum limit. 

This result is consistent with the experimen- 
tal value: ao/mTT = — 1.43(61)(l/GeV^) (aorn,r = 
—0.028(12)). The difference from the predic- 
tion of the current algebra given by ao/rriTr = 
-2.3(l/GeV2) (aom^ = -0.045) is about 1.5cr. 
Since scaling violation is not small, and our data 
points are far from the continuum limit, as seen in 
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Figure 4. ao/m^ at the chiral limit at each (3. 

Fig. U, further calculations nearer to the contin- 
uum limit is desirable. In addition studies with 
the Staggered fermions should be repeated in a 
systematic manner for comparison with present 
results. 
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